Redshift Evolution of the Merger Fraction of Galaxies in CDM Cosmologies 



S. Khochfar and A. Burkert 

Max-Planck-Institut fiir Astronomie, Konigstuhl 17, 
D-69117 Heidelberg, 
Germany 

khochf arOmpia-hd . mpg . de 

burkert Smpia-hd . mpg . de 

ABSTRACT 

We use semi-analytical modelling of galaxy formation to study the redshift evolution 
of the galaxy merger fractions and merger rates in ACDM and quintessence (QCDM) 
cosmologies, their dependence on physical parameters such as the environment, the 
merger timescale, the way major mergers are defined, and the minimum mass of objects 
taken into account. We find that for a given final halo mass the redshift dependence 
of the merger fraction Fmg and the resulting merger rate can be fitted well by a power 
law for redshifts z <1. The normalization -Fmg(O) and the slope m depend on the final 
halo mass. For a given merger timescale tmerg and an assumed maximum mass ratio 
Rmajor for major mergers, Fmg{0) and m depend exponentially on each other. The 
slope m depends logarithmically on the ratio of the final halo mass and the minimum 
halo mass taken into account. In addition, the local normalization Fmg{0) increases for 
larger Rmajor while m decreases. We compare the predicted merger fraction with recent 
observations and find that the model cannot reproduce both the merger index and the 
normalization at the same time. In general the model underestimates Fmg{0) and m by 
a factor of 2. 

Subject headings: dark matter — galaxies: formation — galaxies: interactions — 
methods: numerical 



1. INTRODUCTION 

Estimating the frequency of mergers in the universe is a challenging task. Besides the problems 
of defining a merger in contrast to an accretion event and finding such events, there is also the 
problem of the dependence on the environment, and the estimate of the merger timescales. For 
example, van Dokkum et al. (1999) find that the merger rate evolves as Rmg oc (1 + z)™" with 
m = 6 lb 2 in rich cluster, whereas the merger rate of field galaxies is found to evolve less strongly. 
In a recent study Le Fevre et al. (2000) find m = 3.4 it 0.6 using visually classified mergers and 
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m = 2.7 lb 0.6 using close galaxy pairs in a population of field galaxies. Previous studies found 
m = 3.4±1 (Carlberg, Pritchet, & Infante 1994), m = 4±1.5 (Yee & Ellingson 1995), m = 2.8±0.9 
(Patton et al. 1997), m = 2.01 ± 0.52 (Roche & Bales 1999) and m = 2.1 ± 0.5 (Consclicc 2000). 
This spread in the values of the merger index m is partly due to different methods used in deducing 
the merger rates (see e.g. Patton et al. 1997; Abraham 1999). Correcting for selection effects in 
close pair studies, Patton et al. (2000) estimate that approximately 1.1% of all nearby galaxies with 
— 21 < Mb < — 18 are undergoing mergers. 

On the theoretical side, Gottlober, Klypin, &; Kravtsov (2001) used N-body simulations and 
merger trees based on the Press-Schechter formalism, to derive the merger rate. They found m = 3 
for dark matter halos. In earlier studies of merger rates in N-body simulations Kolatt et al. (1999) 
found m = 3 and Governato et al. (1999) found m = 3.1ib0.2 in a critical universe and m = 2.5ib0.4 
in an open universe. In a previous semi-analytical approach Lacey &: Cole (1993) calculated the 
accretion rate of baryonic cores. They assumed that each halo has only one baryonic core, neglecting 
the effect of multiple cores in a halo. 

The purpose of this letter is a detailed investigation of the galaxy merger fraction and rate. 
In the following section the redshift dependence of the merger fraction and its dependence on the 
cosmological models, on the environment represented by the final dark halo mass, on the merger 
timescale, on the minimum mass of observed objects that would be identified as merger components, 
and on the definition of major mergers are investigated. Besides allowing a better understanding of 
how the merger rates of different observed samples are related, these estimates will test cosmological 
models. 

2. THE MODEL 

We study spatially fiat CDM cosmologies with = 0.3, Qa = 0.7 and Hq = 65 km s~^Mpc~^. 
We also study a quintessence model with Qq = 0.7 and an equation-of-state w = p/pQ = —2/3 (e.g. 
Caldwell, Dave, &; Steinhardt 1998). Merger trees of dark matter halos were generated using the 
method described by Somerville & Kolatt (1999), which is based on the extended Press-Schechter 
formalism (Bond et al. 1991; Bower 1991). This method was modified as proposed by Gottlober et 
al. (2001). The power spectrum is obtained from the fitting formula of Bardecn et al. (1986) and 
normalized by as- We use the expressions derived by Wang &; Steinhardt (1998) for the value of 
(Tg and the linear growth factor. The history of a dark matter halo is followed back in time until 
the masses of all its progenitors fall below a minimum mass of Mmin = lO^^M©. A progenitor with 
mass below Mmin is assumed to inhabit a small galaxy which has 1/10 the mass of the surrounding 
dark matter halo. Whenever two halos merge the galaxies inside of them merge on a dynamical 
friction timescale as described by Kauffmann et al. (1999). 
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3. MERGER FRACTIONS AND RATES 

From the observational point of view one can cither estimate the fraction of visually confirmed 
mergers (Le Fevre et al. 2000) or the fraction of galaxies in close pairs (e.g Fatten et al. 2000, and 
references therein). To deduce the merger fraction it is necessary to correct the observed close pair 
fraction for background/foreground contaminations and to estimate how many of these physical 
close pairs are likely to merge (e.g. Yce &: Ellingson 1995; Le Fevre et al. 2000). Usually one refers 
to the merger rate. The connection between the merger rate Rmg{z) and the merger fraction is 

Rmg{z) = (1) 

where Fjng{z) denotes the fraction of galaxies at redshift z in close pairs which will merge on a 
timescale shorter than tmerg- Since tmerg depends on the separation of pairs, specifying tmerg also 
determines the close pairs. In general, observers measure the separation between galaxies in pairs 
and use the dynamical friction estimate to deduce a merger timescale. We calculate the merger 
fraction by counting the number of galaxies at each redshift which are experiencing a merger on 
a timescale less then tmerg and normalizing them to the total number of galaxies at this redshift. 
The merger fraction of galaxies at redshifts z < 1 is usually approximated by a power law of the 
form: 

Fmg=Fmg{m + zr, (2) 

where F^g(O) is the normalization to the local merger fraction (e.g. Le Fevre et al. 2000). 

For our analysis we consider only binary major mergers, which wc define as mergers with mass 
ratio between Rmajor and 1. Fig 1 shows the result of a representative simulation for a halo of 
mass Mq = 5 X lO^^M© at z = 0, adopting Mmin = 10^*^ M©, Rmajor = 4, and a merger timescale 
of 1 Gyr for the ACDM model. We find in all investigated cases that the merger rate and the 
merger fraction as a function of redshift can be approximated by a power law at redshifts z <1, 
in agreement with the observations. At higher redshifts the merger rate flattens, which was also 
found by Conselice (2000) and Gottlober et al. (2001). 

In general, a range of final halo masses will contribute to the merger events seen in observational 
surveys. To take this into account and to estimate environmental effects wc choose six different 
halo masses Mq at redshift z = (Mq = 10l^ 5 x 10i\ 10^^^ 2.5 x 10^2^5 x in units of 

Mq). Fig. 2 shows the dependence of F„,g{{)) and m, on Mq and tmerg- For increasing Mq, FmgiO) 
decreases and m increases systematically. This trend is consistent with the findings of van Dokkum 
et al. (1999). Varying tmerg corresponds to diflFerent definitions of close pairs. The three curves 
in fig. 2 are exponential laws of the form Fmg{0) = ciexp(c2m), fitted to the merger fractions 
for different tmerg- The parameters used to fit the data points are ci = 0.058 and C2 = —1.23, 
ci = 0.107 and C2 = —1.34, and ci = 0.137 and C2 = —1.42 for tmerg equal to 0.5 Gyr, 1 Gyr 
and 1.5 Gyr, respectively. In the same environment, that is the same final halo mass Mq, Fmg{0) 
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increases with increasing merger timescale as binaries with larger separations are included. The 
merger index m shows only weak variation. 

For computational reasons mergers are only resolved above a minimum mass M^j„. Mergers 
below this mass are neglected. This corresponds to observations with a magnitude limited sample 
of galaxies. The graphs in fig. 3a show the dependence of the merger index m on qm = Mo/Mmin- 
The filled circles are the results of merger trees with constant Mmin = lO^'^M© and varying Mq. We 
compare these results with simulations for constant Mq = 10^^ and varying Mmin, represented 
by open circles. The value of m depends only on the ratio qm SiS m = 0.691n(gM) — 1-77. 

Another important question is the influence of the definition of major mergers on the merger 
rate. The graphs in fig. 3b show the dependence of F^p(O) on different values of Rmajor- An event 
is called major merger if the mass ratio of the merging galaxies is below Rmajor Sind larger than 1. 
As Rmajor increases, Fmg{0) increases. We also find that the merger index m stays roughly constant 
for low Rmajor and decreases at larger Rmajor- A decrease in m with larger Rmajor has also been 
reported by Gottlober et al. (2001). It is a result of the adopted minimum mass for merger events. 
The detectable amount of major mergers with large mass ratios decreases faster with redshift than 
for equal mass mergers, since the small masses drop faster below the minimum mass. In observed 
samples of close pairs Roche & Eales (1999) and Patton et al. (2000) found that Fmg{0) increased 
when they allow for larger Rmajor , which agrees with our predictions. 

How do the theoretical models compare to the observations? The star in fig. 2 is the measured 
merger fraction for field galaxies by Le Fevre et al. (2000), who used Rmajor = 4 and the local 
merger fraction of Patton et al. (1997). They identified close pairs as those which merge on a 
timescale less then tmerg = 1 Gyr. To compare this merger fraction with our estimates one needs 
to take into account that the dark halos of field galaxies can vary over a range of masses and that 
the merger timescale is subject to large uncertainties. We therfore weighted the different merger 
fractions of our sample of field galaxies with halo masses Mq between 5 x IO^^Mq and 5 x IO^^Mq 
according to the Press-Schechter predictions. The merger index m and the local merger fraction 
^mg(O) for different Mq were calculated using the fitting formulae as shown by the graphs in fig. 
2 and fig. 3a. We varied the the range of halo masses contributing to the sample by changing the 
lower bound of halo masses from 5 x 10^^ Mq to 2.5 x IO^^Mq and changed tmerg within the range of 
0.5 - 1.5 Gyr. The results of this reasonable parameter range lie inside the shaded region in fig. 2. 
Results for larger tmerg correspond to the upper part of the region and those for larger halo masses 
lie in the right part of the region. A comparison of our results with the observations shows, that 
the predicted merger index m and the normalization Fmg{0) are a factor 2 smaller than observed. 

As a possible solution to this problem we have studied a quintessence model with w = —2/3. 
The QCDM model shows a shallower increase in the comoving number density of mergers with 
redshift than the ACDM model. There is however not a significant difference in the merger fractions. 
This results from the fact that the difference in the comoving number density is compensated by 
the length of the redshift range contributing to the merger fractions. The QCDM universe with an 
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age of ~ 1.36 x 10 years is younger than the adopted ACDM universe with an age of ^ 1.45 x 10 
years, which is the reason why the same tuierg refers to a larger redshift range in the QCDM universe. 
This result also emphasizes, that it will not be possible to break the degeneracy of these models 
by measuring merger rates. Comparing Le Fevre et al. (2000) results with those obtained from the 
combined Caltech Faint Galaxy Redshift Survey (CFGRS) and Canadian Network for Observational 
Cosmology field galaxy survey (CN0C2) (R. Carlberg, private communication), which includes also 
minor majors, reveals that including minor mergers leads to a smaller merger index m which is 
consistent with the predictions of our simulations. It is therefore not surprising that these two 
observed merger indices differ. 



4. DISCUSSION AND CONCLUSIONS 

Using semi-analytical modelling we recover a power law for the evolution of the merger rates 
and fractions at z < 1, as has been reported in earlier work. Varying the final mass Mq, the 
local merger fraction F^g(0) shows an exponential dependence on the merger index m of the 
form Ffng = ci exp(c2m). The actual values of the parameters Ci depend mainly on the merging 
timescale and on the definition of major mergers. Our predictions that m will increase and Fmg{0) 
will decreases in more massive environments is in qualitative agreement with observations. The 
merger index m depends on the environment through the mass ratio qm- The logarithmic function 
m = C4ln{qM) + C5 fits the data well. We find a similiar behavior as Patton et al. (2000), which 
have shown that if they extend their galaxy sample to fainter magnitudes the local merger fraction 
rises. In addition, we also find that the merger index decreases with higher mass ratios. This is also 
bccing found by comparing the results of the combined CFGRS and CN0C2 sample with those of 
Le Fevre et al. (2000). The adopted QCDM model does not show any significant difference to the 
ACDM model. Therefore it is not possible to distinguish between these two models by measuring 
the merger rate of galaxies. 

Our model predicts values for Fmg{0) and m which are too small by a factor of 2 compared 
with the predictions by Le Fevre et al. (2000) who used the local merger fraction estimate of Patton 
et al. (1997) which was derived with a different definition of major mergers than theirs. As we have 
shown, the definition of a major merger is crucial for the expected merger fraction. Our results 
indicate that the local merger fraction Fmg{0) for the galaxy sample of Le Fevre et al. (2000) who 
used Rmajor = 4 must be less than the value measured by Patton et al. (1997) who used a larger 
Rmajor- A Smaller value of Fmg{0) would lead to an even larger discrepancy in m compared to our 
results. Another issue might be observational errors, like projection effects or uncertanties in the 
merger timescale estimates. Our results strongly emphasize that the comparison of merger fractions 
deduced from different samples and with alternative techniques is questionable if the adopted mass 
range and the definitions of close pairs and major mergers are not the same. 



We thank Ray Carlberg for pointing out the results from the CFGRS and CN0C2 data and 



-6- 



Stefan Gottlober for usefull comments. 

REFERENCES 

Abraham, R. G. 1999, in lAU Symp. 186, Galaxy Interactions at Low and High Redshift, ed. J. E. 
Barnes, & D. B. Sanders (Dordrecht: lAU), 11 

Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 

Bond, J. R., Cole, S., Efstathiou, G., & Kaiser, N. 1991, ApJ, 379, 440 

Bower, R. G. 1991, MNRAS, 248, 332 

Caldwell, R. R., Dave, R., & Steinhardt, R J. 1998, Physical Review Letters, 80, 1582 

Carlberg, R. G., Pritchet, C. J., & Infante, L. 1994, ApJ, 435, 540 

Conselice, C. J. 2000, preprint (astro-ph/0012008) 

Gottlober, S., Klypin, A., & Kravtsov, A. V. 2001, ApJ, 546, 223 

Governato, F., Gardner, J. P., Stadel, J., Quinn, T., & Lake, G. 1999, AJ, 117, 1651 

Kauffmann, G., Colberg, J. M., Diaferio, A., & White, S. D. M. 1999, MNRAS, 303, 188 

Kolatt, T. S. ct al. 1999, ApJ, 523, L109 

Lacey, C. & Cole, S. 1993, MNRAS, 262, 627 

Le Fevre, O. et al. 2000, MNRAS, 311, 565 

Patton, D. R., Pritchet, C. J., Yee, H. K. C, Ellingson, E., & Carlberg, R. G. 1997, ApJ, 475, 29 

Patton, D. R., Carlberg, R. G., Marzke, R. O., Pritchet, C. J., da Costa, L. N., & Pellegrini, P. S. 
2000, ApJ, 536, 153 

Roche, N. & Eales, S. A. 1999, MNRAS, 307, 703 

Somerville, R. S. & Kolatt, T. S. 1999, MNRAS, 305, 1 

van Dokkum, P. G., Franx, M., Fabricant, D., Kelson, D. D., & lUingworth, G. D. 1999, ApJ, 520, 
L95 

Wang, L. & Steinhardt, P. J. 1998, ApJ, 508, 483 
Yee, H. K. C. k Ellingson, E. 1995, ApJ, 445, 37 



This preprint was prepared with the AAS lAT^^X macros v5.0. 



0.020 



0.015 



0.010 



0.005 



0.000 




1.0 1.2 1.4 1.6 1.8 2.0 

1+z 



Fig. 1. — Power law fit to the merger fraction for Mq = 5 x lO^^M©, R 
and t 



ma] or 



4, Mr, 



0) 



merg 



1 Gyr in the ACDM model. The solid line represents the data and the long dashed 



line the power law fit for 2 < 1. 
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Fig. 2. — Showing the merger fraction F^ng at z = versus the powerlaw slope m for major merger 
events with mass ratios less than Rmajor = 4. The data points correspond to different values of 
the merger timescale tmerg and final halo mass Mq. In the upper part of the figure the final halo 
masses in units of Mq are indicated. Halos of the same mass Mq have roughly the same value of 
m. The curves show exponential laws, fitted to the data for tmerg = 0.5 Gyr, 1 Gyr and 1.5 Gyr 
respectively. The shaded region represents the Press-Schechter weighted average merger fraction of 
galaxies in dark halos for the same range of tmerg as mentioned above. The star indicates F^g(O) 
and m as estimated by Le Fevre et al. (2000). The triangle is the result from the combined CFGRS 
and CN0C2 data (R. Carlberg, private communication). 




Fig. 3. — Panel (a) shows the dependence of the merger index m on the mass ratio qm- The points 
are fitted by m = C4 In(Mr) + C5 with C4 = 0.69 ± 0.09 and C5 = -1.77 ± 0.47. Panel (b) of this 
figure shows the local merger fraction for cases with Mq = lO^^m©, Mmin = lO^^'m©, and varying 
Rmajor {Rmajor = 3, 4, 6). Larger values of Rmajor show larger values of Fjng{0)- Both results refer 
to Rmajor = 4. The graphs in (a) and (b) refer to the ACDM model and tmerg = 1 Gyr. 



